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Abstract 

Polymer Quantization (PQ) is a background independent quantization scheme that is deployed 
in Loop Quantum Gravity. This framework leads to a new short-distance (discretized) structure 
characterized by a fundamental length. In this paper we use PQ to analyze the problem of a particle 
bouncing on a perfectly reflecting surface under the influence of Earth’s gravitational field, what we 
have called “XTie Polymer Bouncer”. In this scenario, deviations from the usual quantum effects 
are induced by the spatial discreteness, but not by a new short-range gravitational interaction. 
We solve the polymer Schrodinger equation in an analytical fashion, and we evaluate numerically 
the corresponding energy levels. We find that the polymer energy spectrum exhibits a negative 
shift compared to the obtained for the quantum bouncer. The comparison of our results with 
those obtained in the GRANIT experiment leads to an upper bound for the fundamental length 

o 

scale, namely A -C 0.6 A. We find polymer corrections to the probability of transitions between 
levels, induced by small vibrations, together with the probability of spontaneous emission in the 
quadrupole approximation. 
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I. INTRODUCTION 


One of the main challenges in physics today is the search for a Quantum Theory of Grav¬ 
ity (QTG). On the theoretical front, one of the major problems is that the introduction 
of gravity into quantum held theories appears to spoil their renormalizability, and from 
the experimental point of view, the major difficulty is the lack of experimentally accesible 
phenomena that could shed light on a possible route to QTG. This situation gives rise to 
the possibility of existence of other fundamental interactions in nature [Tj, thus providing a 
strong motivation to the study of gravity at short ranges. The Newton’s law has never been 
tested at distances below of 0.1mm, thus opening the possibility of studying extra finite 
range (non-Newtonian) gravitational interactions. In this scenario different phenomenolog¬ 
ical approaches have been considered , for example extra Yukawa-type forces [2]. Theories 
with large extra spatial dimensions have also been studied. A proposal consistent with string 
theory is due to Arkani-Hamed, Dimoupoulos and Dvali (ADD), in which gravity departs 
from Newton’s inverse square law at scales which could be as large as a millimeter [3jj. ADD 
conjectures the existence of two or more additional dimensions in which gravity, but not 
the strong or electroweak forces, might be acting, diluting itself by spreading its lines of 
force into these extra dimensions. Essentially, this would explain the apparent weakness of 
gravity. 

Considering the progress that has been achieved recently in neutron physics experiments, 
it has been proposed that they can reveal deviations from Newton’s law at short distances. 
For example, Nesvizhevsky and co-workers use the best experimental measurements per¬ 
formed with neutron scattering experiments to give constraints for the hypothetical Yukawa- 
type force [2j, and on the other hand, Frank and co-workers propose that slow neutron 
scattering off atomic nuclei with null spin may provide an experimental test for large extra 
spatial dimensional gravity [4], but no experiment in this regard has been performed so 
far. Note that in these works the authors study the quantum effects caused by an extra 
short-range gravitational potential. 

In this paper we adopt a more fundamental point of view. The high energy behaviour of 
quantum fields is intimately connected with the structure of spacetime at short distances. 
For example, Loop Quantum Gravity (LQG) asserts that continuous classical spacetime is 
replaced by quantum spin networks on small scales. Polymer Quantization (PQ) is a back- 
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ground independent quantization scheme that is deployed in LQG. In this programme the 
notion of discreteness is built-in, and consequently the momentum operator p is not realized 
directly as in Schrodinger quantum mechanics, but arises indirectly through the translation 
operator U\ = e -1 ^. PQ has been used to study quantum gravitational corrections to tem¬ 
poral dynamics [3 IB] and to statistical thermodynamics [7] of simple quantum systems. In 
this paper we use PQ to study the problem of a particle bouncing on a perfectly reflecting 
surface under the influence of Earth’s gravitational held (described by Newton’s law), that 
we have called The Polymer Bouncer. In this scenario, deviations from the usual quantum 
effects are induced by the spatial discreteness, but not by a new short-range gravitational 
interaction. 

Recent high-sensitivity experiments with thermal neutrons free falling in the gravitational 
held, called GRANIT, demonstrate that the energy spectrum in the gravitational held’s di¬ 
rection is quantized [ 8 j. This kind of experiments may offer an opportunity to confront 
observation with signatures of quantum gravitational effects. To this end in section |TT| we 
solve the problem of the polymer bouncer in an analytical fashion, and we evaluate numer¬ 
ically the corresponding energy levels. The Polymer-Schrodinger transition is discussed in 


section |III| In section [TV] we predict an upper bound to the length scale of the polymer 
theory by confronting our results with the maximal precision of the GRANIT experiment. 
In sections [V] and VI we find polymer corrections for both, the probability of transitions be¬ 
tween levels induced by small vibrations and the probability of spontaneous emission in the 


quadrupolc approximation, respectively. Finally, our conclusions are given in section VII 


II. THE POLYMER BOUNCER 

We shall consider the problem of a polymer particle of mass m bouncing on a perfectly 
reflecting surface (at z — 0) under the influence of Earth’s gravitational held: g = —ge z . For 
the purposes of this work we will make the following assumptions. Firstly we consider that it 
is reasonable to keep a constant value for g because the small size of the experiments [9], and 
secondly we assume that the Weak Equivalence Principle is valid. To tackle this problem, 
we restrict the dynamics to an equispaced lattice 7 (A) = {An | n G Z + }. The spectrum of 
the position operator, {z n = An}, consists of a countable selection of points from the real 
line, which is analogous to the graph covering 3-manifolds in LQG. Here A is regarded as a 
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fundamental length scale of the theory. 

As in the standard case, the potential for this problem is defined by 




mg\fi , fJL > 0, 

oo , n < 0, 


( 1 ) 


The wave function V-y = V ; ( z n) of a polymer bouncer obeys the stationary polymer 
Schrodinger equation for the vertical motion along vertical axis 


_ Vv*+i - Vy-i + v VVy 


( 2 ) 


where we have defined the following dimensionless quantities 


£ = 


m\ 2 E 

h? 


v = (7 0 /A) 3 , 


(3) 


( s 2 \ 1/3 

with Iq = {tjEtj) the characteristic gravitational length [10]. The appropriate solution 
for this problem is realized through the following boundary conditions: must vanish 

asymptotically as /i —* oo, and i/jq = 0 because of the mirror at ji — {). 

In most cases, second-order linear difference equations with variable coefficients cannot 
be solved in closed form. Fortunately, for the special case of the difference equation ([2]) 
it is possible to obtain an analytic solution in a simple fashion HU- The most convenient 
approach for finding the solution of difference equations with polynomial coefficients is the 
method of the generating function. Let 


+oo 

/(0= £ VvC", (4) 

fi=—oo 

be the generating function for the solution of eq. ([2]). Note that Q is not a generating 
function in the strict sense but a Laurent series, which is uniformly convergent, and thus 
we can differentiate each term. Multiplying eq. (J2| by and summing for fx 6 Z, the 
following differential equation for the generating function is obtained 


d 

d( 


c 2 — 2(l-e)C + l ' 
C 2 


/ (0 = 0 , 


(5) 


for which the solution is 


+oo 

/ (C) = C 2v(1 -‘ V(H) = c- 2 ” (I - e> 4„(2t.) C", 

fi ——OO 


( 6 ) 
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where we have identified the generating function for the Bessel functions of the first kind 
J n . After relabeling the index in eq. ©. a simple comparison with eq. Q shows that the 
decaying solution of eq. (J2]) is 

Vv = N e J„ +2 u (i_ e ) {2v ), (7) 

where N e is the appropriate normalization factor. By using the recurrence relation for Bessel 
functions, J u+ \ (z) + J v -\ {£) = 2"J u (z), one can further check that eq. ([7]) correctly solves 
the time-independent polymer Schrodinger equation ([2]). This generating function approach 
is not suitable for determining the growing solution of eq. (|2]) , but it is well known that the 
Neumann function also satisfies the recursion equation for Bessel functions. It does not 
however, represent a physical solution for this problem. 

The second boundary condition, ijj 0 = 0, implies that 2v must be chosen as a root of the 
Bessel function J 2u (i- e ) (2v) = 0, where 2v (1 — e) G R (not necessarily an integer). Let j\ n ,r] 
be the r-th zero of the Bessel function J n . Then the condition 


2u j[2v(l—e),n] 


( 8 ) 


defines the quantized energy levels of the polymer bouncer as a function of v. Unfortunately 
eq. ([8]) cannot be solved analitically, but it can be studied numerically for arbitrary values 
of v. Using the quantization condition (J8]) , we can compute the normalization factor in a 
simple fashion. One obtains 

OO 

£ W = JVM+mi-.) ( 2 0 4a-«) (2") = 1. (9) 

fi=0 


where Ja ,0 ^ (z) = dJ Q^ is the derivative of the Bessel function with respect to its order. 
The substitution of N e into the polymer wave function eq. (JT]) establishes the normalized 
polymer wave function. Now we proceed with the analysis of the polymer energy spectrum. 

For the sake of simplicity, let us assume that the gravitational length Iq is a multiple of 
the fundamental length A, i.e. Iq = sX with s G Z + . In table p] we present the first ten 
energy levels of the polymer bouncer as a function of s. We point out that these numerical 
estimates are good approximations for s > 1. The case s = 1 exhibits difficulties arising from 
the polymer behaviour, which we discuss later. Our results show that the polymer energy 
levels display a negative shift with respect to the energy levels of the quantum bouncer. 
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Such negative shift in energy can be understood by studying the dynamics of the polymer 
bouncer. The Heisenberg equation of motion for the operator p\ yields 


dpx i ! pX 

-Jf = J m 9 [UPaJ = -mg cos — > 


-mg (10) 

where the inequality is established using the domain of the momentum in the polymer 


case, i.e. p E [—7rh/\,7r\/h}. In the first order of approximation, eq. (10) becomes p\ ~ 
—mg ^1 — j. So, clasically, an additional p 2 dependent force acts on the particle which 
decreases the energy of the system. 

Quantum mechanically, the negative shift is in agreement with the correction term of the 
regularized p\ operator in the polymer theory for pX <C ft, i.e. 

A 2 


Px 


p 2 - 


12 ft 2 


P 4 - 


( 11 ) 


With the assumption that the fundamental length is very small compared with the char¬ 
acteristic gravitational length, A -C Iq, the energy shift can be roughly estimated. In first 
order perturbation theory, the enery shift becomes 


A e n = 


X 4 


n\p 4 \^n) = 


m 2 A 


2 \ 4 


(i>n\ ( E n - mgzf I-0 n ) . 


( 12 ) 


2 #' r,, ' r 6ft 4 

The averages in this expression can be easily computed by using the quantum-mechanical 


wave function eq. (A7). The final result is 

Ae n = 


(13) 


120s 1 ’ 

which is in a close agreement with the results reported in table [Tj Due to the lack of precision 
in numerical calculations for the case s = 1 in eq. (J8|, we have employed the energy shift of 
eq. (13) in table [I] only for this case. The energy spectrum of a quantum bouncer also has 
been studied in the Generalized Uncertainty Principle (GUP) framework |9]. In this case 
the modification induced by a minimal length Z min on the energy spectrum is 


E^ up ) = —mgl 0 a n + a 2 l 2 min a 2 n , 


(14) 


where a 2 is a GUP parameter. Note that both the polymer- and the GUP-corrections, 
depend on the minimum length and the zeros of the Airy’s function quadratically, but they 
feature a sign difference. Such opposite tendency is also present for the energy spectrum for 
a particle in a box, as pointed out in Ref. [?i|. 
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Another important consequence of the commutator in eq. (10) is the modification of 


the position-momentum uncertainty relation. We know that for two operators A and B, 
the identity B^j = 

m- In the polymer framework the position-momentum uncertainty realtion reads 


A, B ) | J /4 holds irrespective of the quantization scheme 


(Ai) 2 (Ap A ) 2 = j 


COS 


pX 

h 


rf 

4 




(15) 


Note that the correction term in the right hand side is negative definite, thus implying that 


the uncertainty decreases due to the presence of A. Equation (15) suggests that polymer 
quantum mechanics at short enough wave lengths (of the order of ~ 2A) can exhibit classical 
behaviour. This result has a strong resemblance to the one found, on a different setup, by 
Magueijo and Smolin in deformed special relativity |!1|- In this framework, ’t Hooft has 
discussed the possibility of a deterministic quantum mechanics at Planck scale, supplemented 
with a dissipation mechanism, giving rise to the standard quantum mechanical behaviour at 


larger scales pa- Equation ( |15| ) has been used to derive a new mass-temperature relation 
for Schwarzschild (micro) black holes ra- 

We conclude this section by comparing the polymer and quantum density profiles for 
different energy levels and different values of s. Hereafter we denote by the polymer 
wave function of the n-th state. The discrete plot corresponds to the polymer result, and 
the solid line corresponds to its quantum-mechanical counterpart. Figures [la| to lc show the 
ground state for s — 10, 5 and 1, respectively. Figures [id] to [l| show the first excited state for 
s = 10, 5 and 1, respectively. Figures [Tg| to fTT] show the tenth excited state for s = 10, 5 and 
2, respectively. We expect a close agreement between the polymer and quantum-mechanical 
density profiles for the case s>l, however in fig. [l] we observe that this occurs already for 
s = 5. When the gravitational length is of the order of the fundamental length (s = 1) the 
polymer distribution departs from its quantum-mechanical counterpart, and the polymer 
effects become important. From fig. [II] we can infer that the most significant polymer effects 
are for high energies and s = 1, as expected. 


7 













\ n 
s \ 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

10 

0.011686 

0.0204258 

0.0275773 

0.033895 

0.0396679 

0.0450452 

0.0501165 

0.0549412 

0.0595607 

0.064006 

9 

0.0144258 

0.025213 

0.0340387 

0.0418346 

0.0489574 

0.0555915 

0.0618477 

0.067799 

0.073497 

0.0789795 

8 

0.0182553 

0.031903 

0.0430672 

0.052927 

0.0619345 

0.0703228 

0.0782324 

0.0857557 

0.0929579 

0.0998871 

7 

0.0238393 

0.0416556 

0.056226 

0.0690913 

0.080842 

0.0917831 

0.102098 

0.111907 

0.121296 

0.130328 

6 

0.0324385 

0.0566692 

0.0764773 

0.0939613 

0.109926 

0.124785 

0.138791 

0.152107 

0.164849 

0.177103 

5 

0.0466892 

0.0815348 

0.110001 

0.135113 

0.15803 

0.17935 

0.199436 

0.218523 

0.23678 

0.254331 

4 

0.0728877 

0.127199 

0.171511 

0.210558 

0.246155 

0.279241 

0.310382 

0.339949 

0.368205 

0.395345 

3 

0.129331 

0.22536 

0.303481 

0.372143 

0.434588 

0.492495 

0.546873 

0.59839 

0.647516 

0.694599 

2 

0.289409 

0.501951 

0.673219 

0.822395 

0.956849 

1.1875 

1.3125 

1.4375 

1.5625 

1.625 

1 

1.1235 

1.90471 

2.50631 

3.00953 

3.44616 

3.83292 

4.18004 

4.49437 

4.78077 

5.04291 


TABLE I: Estimates of the first ten (rescaled) energy levels of the polymer bouncer as a function 


of s. 


III. THE POLYMER-SCHRODINGER TRANSITION 

In physics, if a new theory is considered more general that the former, it must yield 
accurate results not only in the new scale but at the former one as well. In this sense, we 
know that Newtonian mechanics can be recovered from relativistic mechanics in the domain 
of low velocities compared with the speed of light in vacuum. Regarding the quantum- 
classical transition the problem is more subtle, given that the conceptual framework of these 
theories are fundamentally different DUE!. In the problem at hand one expects that if the 
lattice spacing A is taken to be sufficiently small, the polymer formulation should reduce to 
the Schrodinger representation. However, this is a delicate issue because A is regarded as 
a nonzero fundamental length scale of the polymer theory, and it cannot be removed when 
working in the polymer Hilbert space H po i y [26], no matter how small A is [5]. 

The Polymer-Schrodinger transition for this problem is quite simple. Taking A to be of 
the order of the Planck length, and the value of the characteristic length for the falling slow 
neutrons case / 0 = 5.87/im |18j . one obtains v ~ 10 88 . Note that this approximation is valid 
even if we take A several orders of magnitude larger than the Planck lenght. Therefore, 
the asymptotic behaviour of the Bessel functions for large arguments is required in eq. ([7]). 
On the other hand, the lo > A limit also implies that we should take a very large number 
of points between two arbitrary points, i.e. (i>l. Therefore /r + 2v (1 — e n ) > 2v 1, 
and the asymptotic behaviour of the Bessel function for large orders is also required. The 
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(d) First excited state, s = 10, 
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FIG. 1: Plots of the polymer (blue dots) and quantum-mechanical (red lines) density profiles for 

different energy levels and different values of s. 


asymptotic expansion for the Bessel functions for these conditions is well known 
n 1 and n > x the following approximation is valid 


1 2 (x-n) 

J ’* w ~3 y—r~ 


J 


1/3 


{2 (x — n)} 2 ^ 3 

3^/x 


+ J- 1/3 


{2 (x — n)} 2//3 
3a/x 


For 


( 16 ) 
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Now by using the identity of eq. (]A5|) we find 

Jn (•£) 


/ 2 A 1/3 

/ 2 A 1/3 

- Ah 

- (n-x) 

w 

\xj 


and by setting n = fi + 2v (1 — e n ) 1 and x = 2 u > 1 we obtain 




(n) 


N n 

1/3 


Ai 


v 


z 

lo 


E n 


mgl 0 


( 17 ) 


(18) 


where z — Xfi and E n is defined through eq. (J3]) . From this expression we also recover the 
quantum-mechanical energy spectrum as proportional to the zeros of Airy functions. 

On the other hand, the normalization factor can be approximated by using the Euler- 
Maclaurin formula, 


li=0 


n)\ 


N 2 


v 


2/3 


djiAi 2 


pA _ E n 
l 0 mglo 


N 2 

= -j^Ai' 2 
<■0 


E n 


mgl 0 


= 1 , 


(19) 


where we have used that Ai (x) —» 0 as x —> oo and Ai = 0. Substituting this result 

^n (Z) , ( 20 ) 


into eq. (18) we finally obtain 




Va 


where ^ n (z) is the wave function for the quantum bouncer eq. (A7). Note that the presence 
of A -1 / 2 is consistent with the fact that the Lebesgue measure has dimensions of length while 
the discrete measure is dimensionless m- 


IV. COMPARISON WITH THE GRANIT EXPERIMENT 

The experimental physics of slow neutrons has undergone significant evolution in the 
last decades. Recent high-sensitivity experiments, called GRANIT, performed by V. V. 
Nesvizhevsky et al. at the Institute Laue-Langevin , show that cold neutrons moving in a 
gravitational field do not move smoothly but jump from one height to another, as predicted 
by quantum theory [ 8 ]. They used an intense horizontal beam of cold neutrons directed 
slightly upwards and allowing the neutrons to fall onto a horizontal mirror. By placing a 
neutron absorber above the mirror and counting the particles as they moved the absorber 
up and down, they found that neutrons are measured only at certain well-defined heights. 
The experimental average values of the two lowest critical heights (taken from pjTj) are 

hf xp = (12.2 ± 1.8 sys ± 0.7 stat ) /im, (21) 

hf P = (21.6 ± 2.2 sys ± 0.7 stat ) /nn, 
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while the theoretical heights (h n = —a n lo) with m = 1.67 x 10 2 'Kg and g = 9.806m/s 2 gives 


h\ = 13.7/irn 


/i 2 = 24.0/im. 


( 22 ) 


The theoretical values are therefore located within the error bars. As a consequence of the 
good agreement between theory and experiment, this finding could be used for bounding 
deviations from the standard theory due to an eventual new physical mechanism. In the 
problem at hand, the (negative) energy shift must satisfy the constraint 


|A£„| < A£W 


(23) 


with A E n the energy shift given by eq. (13) and the maximal experimental error. 


The substitution of eq. (13) into eq. (23) produces an upper bound for A given by 

60/ 0 


A 2 < 


-A£ ex P. 


mga z r 


(24) 


Let us discuss the possible bounds for A. For the first two states, we know that A = 
0.102peV and A= 0.051peV [2Tj . With this values eq. (24) yields A < 10 _6 m. This 


result tells us that the effects of the spatial discretization would be largely unobservable in 
the GRANIT experiment, even with an improvement of the experimental precision. On the 


other hand eq. (24) suggests that a better upper bound can be established if the gravitational 


field is intensified by some mechanism. Following this idea, Nesvizhevsky and co-workers 
have considered the centrifugal states of neutrons, which is the quantum analog of the so- 
called whispering gallery wave [22] . They consider the scattering of cold neutrons by a 
perfect cylindrical mirror with a radius of a few centimeters. They found that neutrons are 
affected by a huge centrifugal accelerations of the order 10 5 — 10 ‘ g. Most neutrons entering 
at a tangential trajectory are deviated to small angles. However, some neutrons could be 
captured into long-living centrifugal states |23j. In the limiting case we are considering here 
(A -C /o) our results can be applied to the centrifugal quasistationary states, but not in the 
polymer regime (A ~ l 0 ), because of the (spontaneous or) explicit breaking of continuous 
symmetries [21]. The implementation of this centrifugal acceleration into eq. (24) yields 

o 

A < 0.6 A, which in turn is greater than the mean square neutron radius r n ~ 0.8 fm, thus 

o 

implying that actually this limit should be read as A -C 0.6 A. This is the best bound we 
can establish on the fundamental length of the polymer theory, with respect to the best 
measurements performed in the GRANIT experiment. 
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V. TRANSITIONS INDUCED BY SMALL VIBRATIONS 


It is commonly believed that A is in the order of the Planck length ~ 10“ 35 m (no known 
data substantiates this conjecture), which is 10 24 times smaller than the upper bound we 
established in the previous section. Therefore our results are not yet feasible for the energy 
scales of present experiments, but they motivates the search of other mechanisms that could 
provide a better bounds for the fundamental length. 

Although the quantum states of the bouncing neutron are fundamentally stable, various 
physical effects can induce transitions. In this framework, vibrations and waviness of the 
mirror surface have been analyzed in Ref. |TS]. The purpose of this section is to study whether 
the polymer quantization scheme induces corrections to the probability of transition between 
states that could be significant for detecting quantum gravitational effects. 

Let us assume that the mirror vibrates with a time-dependent height described by a 
time-dependent function h (t) with support on the equispaced lattice 7 (A) = {Xn\n E Z + }. 
By performing the transformation 


n‘ = II-h (() , (i) = W (i), 


(25) 


the modified polymer Schrodinger equation becomes 

d h? 

= ~ 2mA 2 ^P +1 — ^V+i) + rngX^ 1 (j)n' + 


mqXh + ih-^r 
y 2A 



) 


(26) 

where h is the time-derivative and U\ the translation operator. We observe that this equation 
consists of the standard polymer Schrodinger equation for the polymer bouncer (|2]) plus an 
effective perturbation potential 


V 


M ■i h 

mgXn + ih— 
2 A 


Ux - U- ; 


(27) 


Note that the first term provides no transition between quantum levels because it is position- 
independent, but clearly the matrix elements of the second term relates polymer wave func¬ 
tions with different quantum numbers, and can lead to transitions. In first order perturbation 
theory, the probability of the corresponding transition after an observation time T is equal 
to 


Pr 


[T)= V 


v nm (t) dt 


V 

• n 


h 


ke~ i{ - En ~ Em)t,h dt 


(28) 
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with V nm = i Trim 1 where the matrix elements 

OO 

T nm = (F - U-x) 4 


m) 


(29) 


fi =0 


encodes the transition between states. In appendix [B] we explicitly evaluate these matrix 
elements in a simple fashion. The result is 


2 (e„ - Em) T nm = U’!" Vi" 1 ’ - V 2 


(30) 


)i =0 


where e n are the polymer energy levels reported in table [T] and v is dehned in eq. ([3]). Now 
we focus on the quantum regime. In the limiting case A Zq this equation becomes 


V 

< n 


. mg 


UJn 




2l 0 co, 


2 

nm 


(31) 


where c o nm = (E n — E m ) /K is the quantum-mechanical angular transition frequency. On the 


other hand, the integral in (28) makes no polymer corrections to the transition probabilities, 
but it depends on the level of vibration noise in the spectrometer. Therefore, the lowest 
order polymer correction to the probability of transition is proportional to s~ 3 . In the 
GRANIT experiment an accelerometer which is sensitive in the frequency range from 0 to 
500Hz, has been used to study the transition between the six lowest gravitational neutron 
quantum states [IS]. Unfortunately, this implies that the correction term in eq. (31) would 


be strongly suppresed in the GRANIT experiment. However, it could in principle be used 
to establish an alternative bound for A. For comparison with experimental results, we write 
the probability of transition per unit time in terms of the vertical acceleration of the mirror 

3 " 


Pn 


(rrigy 1 

V h ) oA 


l + (-l) r 


9 


IqOJ, 


2 

nm 


Sa (pJnm ) j 


where 


S a (u>) = lim — 
a V > T—>oo T 


rT 2 

he iLJt dt 


(32) 


(33) 


is the acceleration power spectrum [18J. The lifetime for the n-th level can be calculated, 


summing the contributions (32) for all transitions, namely r n 1 = Yl m ^nPn=>-m- I n H ie 
problem at hand, the polymer quantum levels lifetime are given by 

tn. 


TV, = 


1 + t n Q n v~ 3 ’ 


(34) 
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where t n is the expected lifetime of the rt-th neutron gravitational quantum state due to 
vibrations of the mirror (t n ~ 10 4 s) [T5J and 

_ g ^ ^ S a {uj nrn ) 


^n — V fc / 7 

\ h / Iq 


£(-!)' 

m^n 




(35) 


As in the previous section, we can use this equation to establish an upper bound for the 
fundamental length as a function of the best precision measurements, i.e. A t n < A tf^ with 
A t n — t„ — t„. We obtain 

(36) 


A 3 < l 3 0 


tlPr 


where Atn X P is the maximal experimental error in the measurement of t n . Now we perform 
a numerical analysis based on the data reported in [18] for the neutron lifetimes of neutron 
gravitational quantum states due to vibrations of the mirrors. First we assume that the 
maximal error bar in the measurements of t\ = 10 5 s (the larger lifetime) is in the same order 
that for the (3 decay time, i.e. Af^ X P ~ Is. The best precision measurement of the /3 decay 
time of the neutrons is about (878.5 ± 0.7 sta t ± 0.3 S2/S ) s. In order to evaluate Qi, we assign 
the average acceleration power spectrum S a (co nm ) ~ 10 _10 m 2 Hz 3 to all possible transitions. 
A numerical evaluation of the sum appearing in eq. (35) gives fli ~ 10 3 Hz. With these 


values, eq. (36) yields again A < lOun. This upper bound is in the same order as that 
obtained in the previous section when confronted with the GRANIT experiment for thermal 
neutrons free falling in the gravitational field. With an improvement of the experimental 
precision in time measurements, a better upper bound can be achieved. 


VI. TRANSITION RATE OF A POLYMER BOUNCER 


The observation of spontaneous decay of an excited state in the GRANIT experiment 
would be of interest, since it would be a Planck-scale physics effect. This decay rate is too 
low, as estimated in Ref. [23]. Since the spectrum of a polymer bouncer is slightly shifted 


respect to the standard energy spectrum (13), we expect the rate of this decay will change 


as a trace of quantum gravitational effects via the short-distance (discretized) structure at 
high energies. Nevertheless, although the polymer correction to the decay rate is expected to 
be extremely low, it is of conceptual and theoretical interest to estimate it. In what follows 
we discuss this issue, and we present the detailed calculations in appendix [C] 
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The quantum-mechanical transition rate for a bouncer to make a transition k —> n, in 
the quadrupole approximation, is 

T k ^ n = ^ M 2 f iQkm ( 37 ) 

where M p i is the Planck mass, ou kn is the angular frequency of transition and Q kn = 
m ('ijjkl z 2 If/’n) is the quadrupole moment of transition [25]. In the polymer framework, 
the first A-dependent contribution is due to the shift in energy, and the second one arises 
from the polymer quadrupole moment. For the two lowest quantum states, we find that the 
probability of spontaneous graviton emission (at lowest order in A) is 

r 2^i — T2-J.1 ^1 + ■ ( 33 ) 

So, there is a theoretical difference in the probability of espontaneous graviton emission due 
to the spatial “grainy” structure. Although this phenomenon is further away to be detected 
in the laboratory (even the quantum-mechanical result), it is interesting from the conceptual 
point of view because the spontaneous decay of an excited state would be a true Planck-scale 
physics effect. 

VII. CONCLUSIONS 

In the present work, we have investigated the problem of a particle subject to the Earth’s 
gravitational field (described by Newton’s law) but assuming that space is quantized, as 
suggested by various candidates to Quantum Theory of Gravity (such as LQG, String theory, 
and noncommutative geometries). In this scenario, deviations from the usual quantum 
effects are induced by the spatial discreteness, but not by a new short-range gravitational 
interaction as in Refs. M- We have called this problem “ The Polymer Bouncer ”, by 
obvious reasons. We solved the polymer Schrodinger equation in an analytical fashion, and 
we evaluated numerically the corresponding energy levels. We showed that the polymer 
energy spectrum is in a close agreement with the obtained by using perturbation theory, 
supporting the hypothesis that A <C lo in the quantum domain. 

The implications of the introduction of a nonzero fundamental length scale in quantum 
theory are quite profound. For example, there is a belief that if quantum gravity effects are 
taken into account, the Heisenberg uncertainty relations should be modified. In polymer 
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quantization the problem is subtle because the notion of discreteness. Unlike GUP theories, 
in PQ there is no deformation of the algebra of the observables, rather the Hilbert space 
is such that the momentum operator is realized only indirectly through the translation 
operator. In section |TT] we found that polymer quantum mechanics at short enough wave 
lengths can exhibit classical behaviour, resembling to those found, on a different setup, by 
Magueijo and Smolin in deformed special relativity. This finding is also in agreement with 
the’t Hooft proposal about the possibility of a deterministic quantum mechanics at Planck 
scale. 

Experiments with thermal neutrons free falling in the gravitational field offer an oppor¬ 
tunity to test observations with signatures of quantum gravitational effects. In this frame, 
we have established an upper bound for the fundamental length of the polymer theory by 
confronting our results with the best measurements performed in the GRANIT experiment, 

o 

namely A -C 0.6 A. It is commonly believed that A is in the order of the Planck length 
~ 10 -35 m (no known data substantiates this conjecture), which is 10 24 times smaller than 
the upper bound we established. Therefore our results are not yet feasible for the energy 
scales of present experiments, but they motivate for the search of other mechanisms that 
could provide a better bound for the fundamental length. In section [V] we studied the poly¬ 
mer corrections to the transitions induced by small vibrations, and we found an upper bound 


of the same order as that obtained in section IV In both cases, we established a relation for 
the upper bound as a function of the maximal experimental error performed in the GRANIT 
experiment. The improvement of the precision in such experimental methods will produce 


better bounds according to eqs. (24) and (36). For the completion of this work, in section VI 


we briefly studied the polymer corrections to the quantum-mechanical transition rate for a 
bouncer to make a transition k —> n. This phenomenon, although is further away to be 
detected in the laboratory (even the quantum-mechanical result), it is interesting from the 
conceptual point of view because the spontaneous decay of an excited state would be a true 
Planck-scale physics effect. 
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Appendix A: The Quantum Bouncer 


In this section we recall briefly the problem of a particle of mass m bouncing on a perfectly 
reflecting surface under the influence of gravity, that is, of a particle in the potential 


V(z) = 


mgz, 

oo, 


z > 0 
z < 0 


(Al) 


where g is the Earth’s gravitational field strength near the surface. The wave function ip (z) 
of a quantum bouncer obeys the stationary Schrodinger equation for the vertical motion 
along z axis 


K 1 0' 2 iP 


+ mgzip = Eip. 


(A2) 


2m dz 2 

The solution must obey the following boundary conditions, ip [z) must vanish asymptotically 
as z —> oo, and ip (z = 0) =0 because of the presence of the mirror at z — 0. Making the 

/ 2 \ 1/3 

change of variables x = f- — ~^y, where Iq = ( 2 A, ) * s the characteristic gravitational 


length, eq. (A2) can be written as 


d 2 ip 
dx 2 


= xip. 


(A3) 


The general solution to eq. (A3) can be written in terms of the Airy functions, Ai{x) 
and Bi (x). Since the function Bi (x) goes to infinity as its arguments grows, it is not an 
acceptable solution for this problem, where z is unbounded from above. Then, the solution 
is of the form ip [x] = NAi (x), where N is an appropriate normalization factor. The second 
boundary condition, ip (z — 0) = NAi establishes the quantized energy levels of 

the stationary states 

E n = -mgl 0 a n , (A4) 

where a n is the n-th zero of the Airy function. The sequence of zeros of the Airy function 
has no simple analytic expression, but fairly good approximations can be obtained. For 
negative arguments, the Airy function is related to the Bessel functions by 


A* {— y ) — [ Ji / 3 (0 + J - 1/3 (0]) 


(AS) 


where £ = | y 3 ^ 2 . For the zeros Ai (a n ) = 0, one obtains the following approximate solution 


a n « - 


37r 
~2 


1 

U ~ 4 


1 2/3 


(A6) 
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It is accurate within 1%, even for n — 1, and it is exact in the semiclassical limit n » 1. 


Substitution of (A6) into (A4) establishes a good approximation for the energy spectrum. 


The wave function for the n-th (non-degenerate) state reads 


4> n {z) = N n Ai ( a n + — ) 9 (z ), 


(AT) 


where 6 (z) is the Heaviside step function, and N n is determined from the normalization 
condition, 

/*oo /»00 

/ \^ n (z) \ 2 dz = N 2 l 0 / Ai 2 (y) dy = N 2 l 0 Ai' 2 (a n ) = 1. (A8) 


Appendix B: Matrix Elements T m 


To calculate the matrix element T m n of eq. (29), we start with the polymer Schrodinger 


equation ([2]) expressed in terms of the translation operator, 

= (2 - U x - U . A ) V'*” 1 + 


(Bl) 


One can multiply (Bl) by [U\ — U_\J multiply the polymer Schrodinger equation for 
the wave function by (lJ\ — U_\J yip 1 - 1 to sum them and to calculate the sum for /1 e Z + . 

The result is 

OO OO 

2 E O m) 1 = ^ - U—x \ ^ (2 - t/ A - H_ a ) ^ m) + Mm)|(B 


fi=0 


fi =0 


V (c?i - + (n ++ m) 

/i =0 


The left hand side (LHS) can be computed by using the following simple result 

OO OO 

£dr ) tW;r > = 

/ i =0 /2=0 


(B3) 


So we obtain 


OO 

LHS = 2 (E n -E m )J2 d”’ (b - C/-a) V>i m) = 2 ( £ „ - E m ) T„ m , 


(B4) 


H=0 


where T nm is the required matrix element (29). The first term RHS1 in the right hand 


side of eq.(B2) can be studied as follows. We know that to the lowest order in A we have 
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U\ — U_x — — 2Aj(; and 2 — U\ — U_\ ~ —A 2 ^. Therefore the following approximation is 


valid 


RHS1 = J2(Ux~ U- a) Vi n) (2 - Ux ~ U- a) < m) + (n O m), 

[i=0 

1 


(B5) 


^ -j (%-[/_, 




= ~ 4”’ (^A - U. X ) Vo ■ 

Note that in this case RHS1 — — sVi Vi™* because the boundary condition Vo = 0. 
Nevertheless, in the limiting case A <C Iq becomes RHS 1 = —2A 3, ijj' n (0) Tp' m (0), where is 


(m) 


the quantum mechanical wave function eq. (A 2 ). 


For computing the second term RffS2 in the right hand side of eq.(B2) we use the formula 

OO OO 

E /‘d n Aa y ( r ) = E d m) o* t i) & t av'“>. (b6) 

/d=0 /d=0 

The term RffS2 can be written as 

OO 

‘ ~ ~ (B7) 


RHS2 = -u-'Ed"’ (V + & ,) V, 
[1=0 


(m) 

^ ‘ 


After multiplying the Schrodinger equation eq. (B1) by the wave function and perform¬ 
ing the sum for /i 6 Z + , one obtains 

OO OO 

E d n) ( & a + &-a) v<. m) = u- 1 e /‘y ( rV' m) . <bs) 

[i=0 [i=0 

where the orthogonality of the polymer wave functions have been used. Then RffS2 becomes 


iSffS2=-u- 2 EM”V. 

[1=0 


(m) 

M ‘ 


(B9) 


The substitution of (B4), (B5) and (B9) into (B2) establishes eq. (30). 


In the quantum-mechanical regime A <C Z 0 , we approximate the sum appearing in eq. (30) 
by using the Euler-Maclaurin formula, 

™ i r oo 


E 

[i=0 


o Jo 


2 (—i) n ~ m 

Z^n (z) l/) m (*) dz = — - -2 = 2 (-1)' 

yCLn ) 


(Z^) , (BIO) 
V nm J 


where a n is the n-the zero of the Airy function and co nm = (E n — E m ) /h is the quantum- 


mechanical angular transition frequency. As discussed in this section eq. (B4) becomes 
RHS 1 = —2A 3 ^ (0) 7p' m (0) for A lo■ These results establish the required matrix elements 


Rum (eq. (31)) at quantum level. 
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Appendix C: Transition Rate of a Polymer Bouncer 


In this section we present the calculations of the transition rate of a Polymer Bouncer. 
The polymer transition frequency (in the first order approximation), u£ n , is given by 

Ctk T &n 


^kn ^kn ( 1 


60s 2 


(Cl) 


This formula is obtained with the help of eq. (13). On the other hand, we must evaluate 
the polymer quadrupole moment Q% n = rriX 2 To this end we use first order 

perturbation theory. The perturbation shifts the wave function to 


l'MA)) = hM- 


A 2 


24 mh 2 


E 

k^n 


(^klP^ | tpn) 

E n — Eh 


I il>k ), 


(C2) 


where p is the standard quantum-mechanical momentum operator. By using the Hamilto- 
2 

nian, H = ^ + mgz, the first order shift in the wave function becomes 

3 


l*n (A)) = |^n) - (f) E (- 1 )'" 71 F ‘ l^> ’ 

W ljt n 


(C3) 


where we have defined 


Fi n = 


6 


3(Ofc ®n) . (®fc ®n) 

For the evaluation of the expectation value appearing in eq. 
know results 


(C4) 

we have used the well 


\n—k 


2 ( — 1 ) 

(V’fcI Z |t/7i) T 2 t 0 

(®fc ®n) 


In 


(<l>k\ Z 2 \<l>n) = 


24 (-1) 


k—n—1 


\4 *0 


j' 2 

to. 


(C5) 


iflk ®n) 

With the help of eq. ( |C3[ ), the polymer quadrupole moment can be written in terms of the 
quantum-mechanical quandrupole moment {Qkn — m k \ z 2 \4> n )) as follow 

2 


Qkn — Qkn ( ^ 


^ (- 1 )'"* FfkQkn + F ^Q 

. l^k l^n 


kl 


(C6) 


The substitution of eqs. (Cl) and (C6 ) into eq. (37) produces a general expression for the 
transition rate of a polymer bouncer. For the transition between the two lowest quantum 


states, 2 —> 1, we can evaluate numerically the sums involved in eq. (C6). The approximate 


final result is eq. (38). 


20 



















Acknowledgements 


A. Martin-Ruiz would like to thank Valery Nesvizhevsky for the warm hospitality during 
his visit to the Institute Laue-Langevin. AF acknowledges support from DGAPA-UNAM 
and CONACYT. LFU and AMR acknowledge support from the project DGAPA-UNAM- 
IN109013. 


[1] W.-M. Yao et al., J. Phys. G 33, 1 (2006). 

[2] V. V. Nesvizhevsky, G. Pignol and K. V. Protasov, Phys. Rev. D 77 , 034020 (2008). 

[3] N. Arkani-Hamed, S. Dimopoulos, G. Dvali, Phys. Lett. B 429, 263 (1998). 

[4] A. Frank, Piet Van Isacker, Joaquin Gomez-Camacho, Physics Letters B 582 , 15 (2004). 

[5] A. Martin-Ruiz, Phys. Rev. D 90, 125027 (2014). 

[6] A. Martin-Ruiz et al., Rev. Mex. Fis. 61 , 182 (2015). 

[7] G. Chacon-Acosta, E. Manrique, L. Dagdug and H. A. Morales-Tecotl, SIGMA 7, 100 (2011). 

[8] V. Nesvizhevsky et al., Nature 415 , 297 (2002). 

[9] F. Brau and F. Buisseret, Phys. Rev. D 74 , 036002 (2006). 

[10] R. L. Gibbs, Am. J. Phys. 43, 25 (1975). 

[11] M. Ehrhardta and R. E. Mickens, Journal of Computational and Applied Mathematics 172 , 
183 (2004). 

[12] G. M. Hossain, V. Husain and S. S. Seahra, Class. Quantum Grav. 27 , 165013 (2010). 

[13] J. Magueijo and L. Smolin, Phys. Rev. D 67, 044017 (2003). 

[14] G. ’t Hooft, Class. Quantum Grav. 16 , 3263 (1999). 

[15] P. Jizba, H. Kleinert and F. Scardigli, J. Phys.: Conf. Ser. 306, 012026 (2011). 

[16] A. Martin-Ruiz, J. Bernal and A. Carbajal-Donhnguez, Journal of Modern Physics 5, 44 
(2014). 

[17] J. Bernal, A. Martin-Ruiz and J. C. Garcia-Melgarejo, Journal of Modern Physics 4, 108 
(2013). 

[18] C. Codau, V. V. Nesvizhevsky, M. Fertl, G. Pignol, K. V. Protasov, Nuclear Instruments 
and Methods in Physics Reasearch Section A: Accelerators, Spectrometers, Detectors and 
Associated Equipment 677 , 10 (1012). 


21 



[19] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, Edited by A. 
Jeffrey and D. Zwillinger, 4th Edition (Academic Press, New York, 1994). 

[20] E. Flores-Gonzalez, H. A. Morales-Tecotl and J. D. Reyes, Annals of Physics 336, 394 (2013). 

[21] V. V. Nesvizhevsky et al, Eur. Phys. J. C 40, 479 (2005). 

[22] V. V. Nesvizhevsky et al., Phys. Rev. A 78, 033616 (2008). 

[23] V. V. Nesvizhevsky et al, Nature 6, 114 (2010). 

[24] G. Date and N. Kajuri, Class. Quantum Grav. 30, 075010 (2013). 

[25] G. Pignol, K. Protasov and V. V. Nesvizhevsky, Class. Quantum Grav. 24, 2439 (2007). 

[26] The kinemathical Hilbert space can be written as H po i y = L 2 (!Q, d^)> with d/i^ corresponding 
Haar measure, and the real line endowed with a discrete topology 


22 



